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CHAPTER 1 


INTRODUCTION 

The original objectives of this project were to build a test system for the NASA 
255/223 Reed/Solomon encoding/decoding chip set and circuit board. This test system 
was then to be interfaced with a convolutional system at Marshall to examine the 
performance of the concantinated codes. After considerable work, it was discovered 
that the convolutional system could not function as needed. This report documents 
the design, construction, and testing of the test apparatus for the R/S chip set. The 
approach taken was to verify the error correcting behavior of the chip set by injecting 
known error patterns onto data and observing the results. Error sequences were 
generated using pseudo-random number generator programs, with Poisson time 
distribution between errors, and Gaussian burst lengths. Sample means, variances, 
and number of un-correctable errors were calculated for each data set before testing. 

The NASA 255/223 Reed Solomon encoder/decoder is a set of five Application 
Specific Integrated Circuits (ASICs) designed by Advanced Hardware Architectures 
and NASA’s Microelectronics Research Center at the University of Idaho. An 
evaluation circuit board was provided to mount the chips for testing. The encoder 
is a single chip, with three control inputs and the eight data inputs. At the output 
there is a single control line and eight data lines. The decoder consists of the 
remaining 4 chips, using four control inputs and eight data lines. There are two 
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condition flags at the decoder output, along with the data lines. The evaluation 
board has two switches, which provide additional control signals. This report is 
organized into five sections. The first section describes the procedure required to 
operate this chip 6et. Implementing these steps and verifying the result is the 
objective of the test system. 

The second section details many of the operating characteristics of the chip set 
and circuit board. This information was extracted from NASA’s documentation for 
the Encoder Chip, Decoder Chip Set, and Decoder Board. Also helpful were 2 
updated versions of the AHA system documentation and conversations with Dr. Maki 
and Don Wiseman of the University of Idaho. 

The third section covers the peripheral hardware used to support the AHA 
circuit board. To provide flexibility a PC-AT was chosen as data source and system 
controller. Two parallel ports are providing the interface between computer and test 
apparatus. 

The fourth section deals with software, describing the programs for generating 
data, analysis, and control. The procedure to modify this code is also described. 

The final section contains sample test data (with plots), results for the K/S 
system, and program listings. These programs are available on disk, both compiled 
(*.exe) and in ASCII. 



CHAPTER 2 


OPERATING PROCEDURE 

The following is the sequence of events necessary to use the chip set on the 
evaluation board. 

2.A Encoder 

The test program ENCODE.EXE exercises the encoder. 

1 Power on, input off 

2 Clock high 

3 Select Interleaving depth, wait 1.5 microseconds. 

4 Reset Encoder - Bring INC low, clock low/high 32 times, leave clock 
high and INC high. 

5 Input first block of data (between 196 and 223 bytes). The length of 
this block sets the length of all subsequent blocks until a reset. 

6 Stop data, bring INC low, and cycle clock low/high 32 times, leave 
clock high and INC high. 

7 Goto 4 


Data appears at output after a delay of 2 clock cycles. OUTC is simply 
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INC passed through along with the data. It is high when data is output 

and low when the 32 check bytes are output. It will be necessary to generate a Block 

Begin Not signal to tell the decoder where data blocks begin. 

2.B Decoder 

1 Power on, input off 

2 Clock high 

3 Reset Decoder - Bring SystRstN low for 8 clocks, high for at least 
7 clocks. 

4 Bring BlkBegN low for one clock cycle coincident with the first data 
byte input. 

6 Clock in the rest of the code block. 

6 Goto 4, repeat until buffers are empty. 

Output is valid when DataVldN is low. If ErrPattN is high, the corrected data 
along with a status byte are output. The status byte is in place of the first check 
byte, the remaining 31 check bytes are set to zeros. 

If ErrPattN is low, the detected error pattern is supplied at the data output 


lines. 



CHAPTER 3 


SYSTEM CHARACTERISTICS 

This section details many of the operating characteristics of the chip set and 
circuit board. This information was extracted from NASA’s documentation for the 
Encoder Chip, Decoder Chip Set, and Decoder Board. Also helpful were 2 updated 
versions of the AHA system documentation and conversations with Dr. Maki and Don 
Wiseman of AHA 

NOTE: - Always turn on power before applying an input signal. 

3.A Encoder 

3A.1 Input Description 

- All inputs are TTL compatible 
DATA - TTL input Data byte. 

Clock - Rising edge triggered, up to lOMhz, data must be valid for 33ns 
before the clock rises and held for 10ns afterwards. 

INC - Goes high with first data byte, stays high during data input. 

- The number of clock cycles that this lina is high dete rmin es 
block length. The decoder allows the block length to vary 
between 255 and 228 bytes, inclu ding the 32 check bytes. 

- To c h a n ge block length, you must reset the encoder and decoder. 
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- Don’t forget to clock out the data in the decoder pipeline. 
SELECT - Sets interleave depth, currently set to 1. Lines must be 

stable 1.5 us after "power up sequence" (reset) or change in 
interleave depth. 

T imin g - Clock can stop & restart, as long as the data flow remains 
synchronized with the clock. 

- Clear encoder on start up: bring INC low, clock 32 times. 

- There is a 2 clock delay for data through encoder. 

- At end of each data block, stop data flow and hold INC low for 
32 clocks to output check bytes. 

• The INC signal can be used to trigger an oscilloscope on 
the first data byte or the first check byte. 

3.A.2 Output Description 

- All outputs are TTL compatible 
DATA - Data and check bytes 

OUTC - INC passed through with the data, goes high with first data 
byte out, low during check byte output. Logic is inverted 
from BlkBgnN, and does not provide the needed one clock strobe. 
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3.B Decoder 

3.B.1 Input Description 

- ErrPattN and SysRstN are CMOS only! 

- BlkBegN - Goes low for one clock cycle at the beginning of each 
code block. The number of clock cycles between the first and second 
strobes of BlkBegN after reset determine the length of the code 
blocks. 

- BlkBegN is a mixed CMOS/TTL signal 

Clock - Separate from encoder clock. Same description as for encoder. 

DATA - Data and check bytes 

DUAL - Select dual or normal field (PROM’s contain the field). 

It is selected by the jumper labeled SI. 

ErrPattN - Error Pattern Not **CMOS** Selects whether the corrected data 
or detected error pattern will be output. 

- Low means error pattern and check bytes are output. 

High means data (correct or not) and status byte are output. 

SysRstN - System Reset Not **CMOS** 

Timing - Reset to clear pipeline at startup or when hang in g block 

length. Bring SysRstN low for at least 8 clocks, then high 
for 7 clocks. Any data in the FIFO will be lost. 

- BlkBgnN must go low for one clock cycle in conjunction with the 
first data byte of every block. 
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- There is a pipeline delay of about 3.5 code blocks between data 
entering the decoder and the corrected results co min g out. 

The delay for a 255 byte block is 908 clocks. There does not 
appear to be any good way to calculate the delay. 

3.B.2 Output Description 

- All outputs are TTL compatible. 

Correct - High indicates that the block of data being output was correctable. 

- Low indicates that the data was not correctable, no 
corrections were made. Whatever was received will be passed 
straight through to the output without change. 

- This Correct is valid when DataVldN is low. 

DATA - Data and check bytes 

DataVldN - Data Valid Not is low when data (correct or not) is being output. 
High for status byte and 31 padded blanks. 

Status Byte - First byte output after data block, corresponds to first check 
byte in timing. Bit 7 = 0 means data valid, bit 7 = 1 means 
uncorrec table errors occurred. Bits 0-4 indicate the binary 
count of corrected errors, up to 16, with bit 4 the most 
significant. The remaining check bytes are all zero. If no 
errors occur, the status byte and 31 check bytes will all be 
zero. 

T imin g - The DataVldN signal can be used to trigger an oscilloscope on 
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the first data byte or the first check byte. Due to time 
variations of the program loops, the scope trace tends to 
j um p around a lot. For some test runs all maskable interrupts 
were disabled, so the only way to stop the program is to reset the 
computer. 

3.C Buffers 

The buffers isolate the ASIC chips from the input and output sockets. 

Switch 2 jumpers the enable line on the 74LS244s to +5 volts or ground. Ground 
enables the buffers, +5 disables them. The 74LS244s have three-state output latches, 
but the control lines are not implemented. The buffers add two 18ns delays through 
the circuit. This does not change the required setup and hold times, but if the board 
is inserted into another synchronous circuit, this may slow the maximu m clock speed 
of that circuit. 

Notes: - The 74LS244 decoder output buffer is wired backwards on the circuit board. 

It must not be connected. 

- Text on the block diagram in the circuit board documentation reads 
upside down. Pin one on the buffers and PROMs is to the top of the 


board. 
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3.D FIFO 

The FIFO is IK x 9, id number "IDT 7202SA35j" in a PLCC package. It is 
available from: 

Integrated Device Technology, Inc. 

4930 Corporate Drive, Suite 1 
Huntsville, AL 35805 205-721-0211 

Sales Engineer Steve Gary provided two of them at no charge. These are CMOS 
chips. 



CHAPTER 4 


PERIPHERAL SUPPORT FOR R/S BOARD 

Figure 1 shows a block diagram of the experimental equipment including the 
computer, error injection, and R/S board. Function of the various subassemblies is 
discussed in this section. 

4.A Computer 

A PC-AT type computer was selected to be the controller and data source for 
this project. The machine used operated at 10 MHz, and had 640K bytes of RAM, 
a 20M byte hard disk, EGA monitor, an 80287 math coprocessor, and two parallel 
ports. The search for fast I/O ports turned up several boards, however none were 
designed for the speeds needed for this project. Examination showed that direct port 
control using assembly language and/or C was essential to ob taining maximum speed 
when downloading data. None of the available I/O boards indicated that they could 
operate faster than the built-in parallel ports operated via assembly language. The 
computer ports provide all timing and control signals. The process used to generate 
computer code to exercise the board via parallel ports is given in the computer 
section. A 2 Megabyte RAM disk was added during testing. 
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4.B Power 

Power supply should be between 4.75 and 5.5 volts, and must not exceed >0.3 
to 7.0 volts. Power must be supplied to the ASIC chips before signals are applied. 
If the power is cycled o£E7on, allow all capacitors to discharge and remove all inputs 
before tur ning back on. Failure to do so may lead to latchup and chip failure! 
Current to the board should not exceed 1.2 amperes. 

On this test set, a 7805 voltage regulator steps down the PC’s 12 volt 8 amp. 
(yellow wire) line to 5 volts for the AHA-NASA circuit. It also limits the current to 
no more than 1.6 amps. The purpose of the 7805 is to prevent the accidental 
connection of the 12 volt line (or any other voltage) from the unkeyed power socket 
to the board. If the 7806 is pushed beyond its limits (in excess of 1.5 amps until it 
reaches thermal overload), it can oscillate. Limiting the current to 1 amp or less 
using a fuse and placing a 0.1 uF capacitor across the output prevents the conditions 
leading to oscillation. A switch can cut off the 12 supply to the 7805, removing the 
power to the R/S board. This is intended for use with the buffered input 
configuration only. DO NOT APPLY POWER TO THE BOARD WITH SIGNALS 
ALREADY PRESENT AT THE UNBUFFERED INPUTS. THIS MAY BE FATAL TO 
A CHIP. A 5 J. volt 1 watt zeiner diode is present to trap overvoltage surges. 

The computers 5 volt 20 amp. (red wire) line supplies all of the other circuits 
used in the test bed. A 2 amp fast blow fuse and a cutout switch protect those 
circuits. A voltage monitor circuit has been designed to shut down the system if the 
supply voltage on the board should exceed the specified safe range. The circuit was 



13 


not used with the testbed, since the computer supply appeared to have good 
regulation. 

4.C Parallel Ports 

The parallel ports are generally known as Centronics ports. This is 
misleading, they are actually somebody’s implementation of Epson’s modification to 
Centronics modification to the Intel 8255A Programmable Peripheral Interface. Not 
all of the modes defined for the Interface are implemented. By experimenting with 
defined modes and writing to individual bits, the location and means of accessing the 
pins were obtained. 
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TABLE 4.1 

Epson/Centronics Printer Port Pinout 


pin# 

Standard Definition 

Direction 

1 

Strobe 

Out 

2 

Data 

t» 

3 

tt 

it 

4 

it 

it 

5 


it 

6 


it 

7 

ti 

if 

8 

w 

tt 

9 

it 

n 

10 

AcknowledgeNot 

In 

11 

Busy 

In 

12 

Paper Empty 

In 

13 

+Select 

not defined 

14 

Auto Line Feed 

Out 

15 

Error Not 

In 

16 

Initialize Not 

Out 

17 

-Select Not 

Out 

18-25 

Ground 
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The port is logically divided into four 8 bit ports as follows 

TABLE 4.2 

Parallel Port Addresses 
Port PC/AT Address (Hex) 



Video Card 

Printer Card 

Direction 


(lptl) 

(lpt2) V 


Data 

0X0278 

0X0378 

out/in 

Status 

0X0279 

0X0379 

in 

Print control 

0X027A 

0X037A 

out 

Port control 

0X027B 

0X037B 

— 


The port modes are supposed to control direction and functional relations 
between the logical ports. In particular, in modes 1 or 2, data appearing on the 
output pins should auto matically be accompanied by a strobe on pin one. The speed 
of the port was determined by writing a short assembly routine to dump alter nat i ng 
ones and zeros to the data lines, and measuring the resulting output signal with an 
oscilloscope. Without any hnnHahalriTig overhead, the maximum rate achieved to date 
is about 600,000 bytes per second. This is close to the maximum rate that this 
particular I/O card can physically support, the signals are getting rounded and are 
not reaching their full amplitude. Nothing in the documentation appears to indicate 
that, the port should not be capable of running at the computers bus speed (usually 
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8 Mhz), but your standard cheap I/O card may not be that good. The built-in strobe 
may not be reliable for a high speed clock. The strobe is by far the most rounded of 
all of the output signals. It seems this is intended to delay the strobe until the data 
setup requirements have been met. Simply putting a long RC timeconstant on the 
strobe provides the delay. 

4.D Cables 

The connector on the printer port is a 26 pin D shell. The current cable is a 
36 wire ribbon. Data and ground line alternate, with all of the grounds joined at 
both ends. At higher frequencies it may be desirable to leave one end of all but one 
of the grounds disconnected. The port pins are currently assigned as follows: 

TABLE 4.3 

Ribbon Cable Connections 


pin# 

LPT2 

Function Direction 

LPT 1 

Function Direction 

1 

Not used 

Out 

Not used 

Out 

2-9 

Data 

Out 

Error 

Out 

10 

Not used 

In 

Not used 

Out 

11 

M 

M 

N 

tt 

12 

It 

m 

m 

it 

13 

W 

tt 

« 

it 

14 

Clock Out 

Out 

Clock 

Out 

16 

Not used 

Not defined 

Not used 

Not defined 

16 

Not used 

Out 

Not used 

Out 

17 

INC 

Out 

INC 

Out 
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4.E Error Injection System 

The error signal is generated by the second parallel port on the PC. Data is 
Exclusive ORed with the error signal at the output of the encoder. Note that the 
data has a two clock latency passing through the encoder, so errors are injected onto 
data which entered the system two clock cycles before the error. 

4.F Block Begin Not Generator and Clock Inversion 

The control signal BlkBgnN is generated using the NAND of INC (delayed one 
clock) and OUTC. The signal is delayed one more clock, then applied to the decoder 
input. The timing and logic diagrams for generating BlkBgnN are illustrated in 
Figure 4.1. The clock for the Decoder is also inverted at this stage. This is to 
provide a delay (for setup time) between the application of the clock to the encoder 
and the decoder. This is known as a split phase clock. 

4.G Miscellaneous Notes 

All chips mount with the notch to the top of the board. The Block diagram is 
misleading. All of the ground pins are on the outside edge of the connectors. Don’t 
try to use both the buffered and unbuffered inputs at once. 

The encoder and decoder clocks are not common on the board. 

CORRECT and bit 7 of the status word should indicate the same error 
condition. CORRECT is available during the entire data bl ock, but if it indicates an 
uncorrectable error it will be low during the output, hi g h during the check 
bytes. To latch the check byte, pasB the data through a Hnnlc of D flipflops, then 





D Flip Flops 


Figure 4.1 Block Begin Not Generator and Timing Diagram 
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AND the dock with DataVldN. One dock after the check byte is output DataVMN 
goes high, and the status byte will be on the flipflop output lines. The output of the 
AND signals that the check byte is on the D flipflop output te rminal s. The status 
word may be stored along with the data for a record of how well the system was 
performing. The count of corrected errors might be used as a measure of the 
performance improvement produced by other factors in the system such as addition 
of soft dedsions or a convolutional code. The documentatipn does not indicate when 
ErrPattN can change. 

Maximum rated clock speed is 10 Mhz, yielding an equivalent input bit rate 
of (N-32)/(N) * 10 6 * 8 where N is the number of bytes in a code block, 228 <= N <= 
265. A interesting problem occurred due to the fast response of these chips. The 
counter used to monitor uncorrectable errors was highly erratic. Adjusting the 
sensitivity reduced the magnitude of the miscount, but could not stabilize the 
counter. Connecting two counters in parallel did fix the problem. It seems that the 
rise time on the decoder output is so rapid that the counters had difficulty keeping 
up with the signal. Placing two counters in parallel supplied enough load that the 
signal was damped slightly. The two counters subsequently produced accurate and 
consistent results. 

The signals injected into and produced by the chip set were monitored using a 
dual tdmebase oscilloscope. Using the external trigger on the block begin not or data 
valid not signals allows examination of any part of the waveform. High impedance 
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tuned oscilloscope leads were used to minimize loading effects on the circuit during 
testing. 

The doc um entation indicates that the encoder can implement some form of 
interleaving. It does not provide details of how the interleaving is arranged or how 
it is to be undone. Apparently, only the check bytes are affected. The three 
interleaving select lines are tied low through 4700 ohm resistors, selecting an 
interleave depth of one. Further information has been requested several times, but 
has not been forthcoming. 

Figure 4.2 shows the arrangement of the test apparatus for loopback testing of 
the Reed/Solomon board alone. 











CHAPTER 5 


SOFTWARE 

There are several distinct tests to be performed to satisfy the stated objectives 
of this work. For each of these tests, the statistical distribution, burst le ng th, 
spacing, and other characteristics of the noise and resulting error patterns need to 
be specified. This chapter describes test procedures for the 255/223 Reed/Soloxnon 
encoder/decoder made by AHA for NASA Each 255 byte code block (word) is 
currently arranged into 223 consecutive data bytes (8 bits/byte), followed by 32 
consecutive check bytes. Programs were written in the Microsoft Quick C editor, 
compiled to an assembly language listing with Microsoft’s C Compiler 5.1, optimized 
in assembly using a program editor, compiled to a .obj format using Microsoft’s ma c r o 
assembler MASM, and compiled to .exe format using the C compiler again. This 
involved process is only necessary for mwimnm speed on the output routines. The 
compiling steps have been largely automated using batch files. Ordinary C or 
Fortran code may be used to generate the data and place it into arrays or files. 
These assembly I/O routines can be inserted into C programs using the latest version 
of Microsoft QuickC. Some of the advantages of this method are: 

-The QC (Quick C) editor checks syntax for you. 

-The C5.1 Compiler supports very large data arrays and takes care of 
operating system bookkeeping. 
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•Some simple 6teps in assembly can vastly increase speed. 

•Assembly language gives direct access to the bits of the ports. 

Error patterns were applied to the all zeros data block and/or alternating 
ones/zeros byte data blocks. Testing was performed in two major sections, 
corresponding to the method used in generating the data and error vectors. Initial 
work was performed using "immediate" data, meaning that the data were part of the 
assembly language programs. This is the simplest and fastest implementation. 
Output was monitored using an oscilloscope triggering on the DatavalidN line. Both 
the data and the Correct lines were monitored. The Status byte was also examined 
and in all cases identified the number of correctable errors. This mode was 
principally used to debug the system. The second mode of testing was to generate 
large data files on disk, then dump the data to the output port. 

The clock rate is typically about 225,000 Hz for the parallel immediate mode data 
(1.8 MHz if converted to serial data, this exceeds the data rate of the OMV system, 
but it is not difficult to slow down). Exact speeds are code dependent, but all were 
close to this rate. This corresponds to over 810 million bytes or 6.48 billion bits per 
hour. 

6.A Immedia te Mode Data Testing 

The first testing was primarily aimed at debugging the computer interface with 
the R/S board. These same error vectors were applied to the board once the interface 
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was functioning properly. The data is part of the assembly language program, loaded 
into the CPU registers directly. 

5.A.1 Single Bit Errors 

Tnitinl testing was performed by inserting a single bit error on bit zero of the 
third from last data byte of each code block (byte 221). The program ran in a 
continuous loop clocking out this error on the all zeros data block. Subsequently this 
error was moved to each bit location of that byte. This error pattern was applied for 
about 40 hours over a period of two weeks during initial testing. Appro xim a t ely 6 
hours of this time the system was closely monitored, the balance was to allow burn-in 
of the test system. No indication of incorrect decoding was ever observed. (Incorrect 
decoding is failure to correct errors, failure to identify uncorrectable errors, or 
insertion of additional errors.) 

5.A.2 Bvte Errors 

Most of the testing has been performed using the case of every bit of every error 
vector being all ones. (ie. all error bytes = 11111111). First, the same byte as above 
was corrupted, then data bytes 219 and 221 were both corrupted. Then this pattern 
(bad byte, good byte, bad byte), was moved to the first data location. This pattern 
was used to examine the alignment of the DatavalidN and BlockbeginN timing 
signals. This pattern was applied for about 20 hours over several days. Again no 
faulty decoding was observed. 
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5.A.3 Block Errors 

Next, blocks of data were corrupted, and the error location in the data block 
moved through the entire data block. This was implemented by using a set of nested 
loops adj usting the loop counters. The error block length ranged from a single 
byte up to 17 consecutive bytes. (All error bytes for burst errors were 11111111, 17 
consecutive error bytes isl7x8 = 136 consecutive ones). These processes ran for 26 
hours. All cases of 16 or fewer errors were corrected. The case of 17 errors was 
properly identified as uncorrectable. Lastly, the extreme case of an entire data block 
was corrupted, clocking out all ones errors on all of the data and check bytes. The 
system correctly indicated an uncorrectable error. 

5.B Data File Testing 

A more sophisticated method is to read the necessary data from a hard disk file 
into RAM and then output the data from RAM to the encoder. Tests can be 

dofinod using various error distributions such as Gaussian, Poission, etc. 

The next step in testing was to generate large error data files and run them 
through the system. Two methods are being used to generate these data files: First 
a Fortran program was used to generate Gaussian distributed random numbers. To 
generate about 14 errors per 255 byte code block, find the area corresponding to 
14/255 = 0.0649 under the standard normal curve. This corresponds to the X axis 
value -1.6, so each Gaussian number is interpreted as an all ones error byte if the 
random n umb er is less than or equal to -1.6, or an all zeros error byte if the number 
is greater than -1.6. The decimal equivalents 255 or 0 were stored in a file to be read 
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and output to the encoder. A file containing 25,500 of these error bytes was created 
and is currently r unning on the system. The number of errors per code block can be 
changed by finding the X avis point corresponding to the area under the standard 
normal curve, then re-evaluating the decimal random n umb ers. 

The statistical analysis software SAS was to have been used to generate data files. 
Unfortunately, data format information for the storage of eight bit unsigned integer 
data was not readily available. A F ortran program generating Gaussian random data 
has been implemented following the appendix on random n umb ers in Orfamdus 
Signal Processing book. The data is random burst errors (all ones) with a Gaussian 
burst length (ie. number of bytes/burst) and Poisson time distribution between bursts. 

The data is analyzed before the tests so that the number and length of correctable, 
uncorrectable, and undetectable errors is known for each data file. Co unting the 
results out of the decoder provides collaboration. These same error files could be 
injected into the serial data stream for the concantinated codes. 

Speed problems arise when using the disk file data files. First, the clock out is 
stopped every time, the disk is accessed for more data. The clock rate while data is 
being transferred to the encoder is essentially unchanged, but when memory is 
empty, more data must be read in from disk The rapid data output rate leads to 
constant disk access. Implementation of a 2 megabyte RAM disk to supplement the 
hard disk improved data access and reduce wear on the hard drive. Data 
compression has also been introduced by storing only the length of error bursts and 
the number of clear bytes between error bursts. For the best case, this represents 
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each 256 byte block in two eight bit numbers, for a compression of over one hundred 
to one. The data are stored in two files, along with a discriptor file conta inin g the 
number of terms to read from each data file per memory update. 



CHAPTER 6 


RESULTS AND RECOMMENDATIONS 

Histograms, sample means, and variances of each data set were generated to 
document the characteristics of test data. The following three tests are 
representative samples of the general types of data used and results obtained. The 
results show that the system is working exactly as predicted. For example, in test 
12, there were 42,519 code blocks (86,738,760 bits) known to be uncorrectable out of 
200,892 code blocks (409,819,680 bits) in the test run. For three successive runs the 
experimental count showed exactly 42,519 uncorrectable code blocks. Observation 
of the correctable errors into and out of the decoder (on an oscilloscope) showed that 
these errors were in fact corrected. This means that the errors con taminating 
157,481 code blocks were in fact identified and corrected. 

In the experimental data some code blocks have multiple error bursts, while 
others have no errors. This causes a difference in the number of error bursts and the 
number of code blocks in the tests. 

The histogram data files for the Gaussian Burst lengths are included in the 
Appendices. Histogram plots of both Gaussian Burst length and Poisson spacing 
data are also included. Data files for Poisson contain several thousand entries, so 
they were not included. 
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6.A Sample Tests 

The fo ll o wing three tests are representative samples of the type of experiments 


performed and corresponding data. 


Test 4 

Gaussian Burst: Mean = 1, Variance = 7 Poisson Spacing: Lambda = 128 

Number of: 

error bursts 100,000 

code blocks 99,996 (26,498,980 bytes, or 203,991,840 bits) 

correctable code blocks 99,281 (26,316,666 bytes, or 202,633,240 bits) 

uncorrectable code blocks 716 (182,326 bytes, or 1,468,600 bits) 

This is 0.71% uncorrectable errors. The test was repeated 3 times. The number of 
measured uncorrectable blocks was 716 for each ran, matching exactly the known 
number of uncorrectable errors. 


Test 9 

Gaussian Burst: Mean ** 8, Variance = 4 Poisson Spacing: L ambd a * 128 

Number of: 

error bursts 200,000 

code blocks 104,284 (26,592,420 bytes, or 212,739,360 bits) 

correctable code blocks 66,740 (16,763,700 bytes, or 134,109,600 bits) 

uncorrectable code blocks 38,664 (308,432 bytes, or 2,467,466 bits) 
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This is 36.9% uncorrectable errors. The test was repeated 3 times. The number of 
measured uncorrectable blocks was 38,654 for each run, mat c hing exactly the known 
number of uncorrectable errors. 

Test 12 

Gaussian Burst: Mean = 10, Variance = 8 Poisson Spacing: Lambda = 256 

Number of: 

error bursts 200,000 

code blocks 200,892 (51,227,460 bytes, or 409,819,680 bits) 

correctable code blocks 158,373 (40,386,115 bytes, or 323,080,920 bits) 

uncorrectable code blocks 42,519 (10,842,345 bytes, or 86,738,760 bits) 

This is 21.1% uncorrectable errors. The test was repeated 3 times. The number of 
measured uncorrectable blocks was 42,519 for each run, matching exactly the known 
number of uncorrectable errors. 

6.B PmhlflmB & Suggestions 

The following recommendations break down into two parts: modifications in 
the chip set itself, and information for the user. The ’problems’ described are things 
which make the encoder/decoder set more difficult to use. While a given user can 
adapt his/her system to deal with these characteristics, the chip set could be easier 


to use. 
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1 Problem: Verifying tbe correct function of tbe encoder/decoder independent of 
t h e communicatio n system in which it is embedded is difficult or impossible. When 
building a c ommunica tions system, it is a good idea to test each subsystem 
independently. It may also be useful to be able to request such a self-test of an 
orbiting system. 

Suggestion: Facilitate the performance of a loop-back self-test in the chip set. 

One way to implement a self-test is to feedback the output of the encoder into 
the input of the decoder. This requires generation of various control signals. The 
following two features would make it easier to implement such a self-test. 

* Tn c bide a Block Begin Not generator on the encoder chip. A simple circuit to 
generate this signal is shown in Figure 4.1. 

* Supply a delayed clock signal at the output of the encoder. This is to allow setup 
Hma for data and control input to the decoder. 

Another way tests of the encoder and decoder might be performed is to provide 
an input sequence to each fr om ROM, then compare the encoded/decoded data with 
known good results also stored in ROM. 

2 Problem: The system in which the chip set is embedded must operate the 
encoder and decoder by providing timing and control signals. 

2,A The user must reset both encoder and decoder manually. The present system 
also gives no indication of when a successful reset has occurred. 
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2.R The system using the chip set must provide the encoder and decoder control 
signals INC and Block Begin Not. Generating control signals to tell the chip set 
when every data block begins and ends demands additional work from the 
communication system. 

2.C The encoder input data stream is not continuous, being interrupted every code 
block for the time required to clock out the check bytes. Halting the data flow into 
the encoder while the check bytes are clocked out requires that the data source count 
out each data block individually. 

Suggestion: Make the code/decoder chip set as independent as possible of the 

surrounding system. This might be accomplished by modeling the interface logic on 
a standard co mmunic ations protocol such as RS-232C. The data source asserts 
Ready-To-Send (RTS), the receiver acknowledges with Clear-To-Send (CTS). When 
thu receiver needs to stop accepting data, it negates CTS. The internal function of 
the encoder/decoder would be essentially transparent to the user. 

2A A simple sequential circuit could execute a system reset in response to a 
request. If each chip provided a reset signal, a logic circuit could generate CTS. 

2 R INC could be replaced with two variable length counters, one for the number 
of data bytes, and the second for the number of check bytes. These count lengths 
could be selected using jumpers for standard code lengths. After a system reset, the 
first clock would correspond the first data byte. 
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The decoder block length should also be set using counters as described above 
for the encoder. Once the data stream is initialized, the encoder would know where 
each new block begins without having to be told. 

2.C Making the input data flow continuous requires that the data be clocked out 
of the encoder faster than it is clocked in. If the input data is buffered in a FIFO, 
the input could continue uninterrupted while the check bytes are output. The 
necessary clock rates are related by the factor Code-Block-Length-Out/Data-Block- 
Length-In, 255/223 (=1.143497 , a continued fraction) for the nominal code. 
Multiplying by 255 and then dividing by 223 produces very high intermediate 
frequencies. At a 10 MHz output rate, 2230 MHz intermediate clock rate would be 
necessary. The input data rate would be 223/255 * 10 MHz, or about 8.745 MHz. 
The integer 223 is prime, so the ratio can not be reduced to smaller multipliers. The 
ratio might be adjusted by adding bytes between code blocks (sync, patterns perhaps), 
by shortening the data blocks, or some combination of the two. Perhaps a single low 
frequency reference could be multiplied by 255 and 223 separately to produce the 
needed clocks. 

At 10 MHz, even a small percent error would quickly lead to FIFO overflow or 
underflow. Allowing one of the clocks to vary its period slightly (within bounds 
determined by the FIFO size) could get around this problem. 

If the above modifications were implemented (except 2.0, the following input 
and output signals would be needed: 
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Encoder 


Input Signals 

Data 

Clock 

RTS (Request-To-Send) 
RR (Request Reset) 


Output Signals 
Data 

Delayed Clock 
OutC 

Block Begin Not 
CTS (Clear-To-Send) 
RV (Reset Valid) 


Decoder 


Data 

Clock 

RTS (Request-To-Send) 
RR (Request Reset) 


Data 

CTS (Clear-To-Send) 
RV (Reset Valid) 
Data Valid Not 
Correct 


3 Provide up-to-date documentation to the end user. A lot of time could have 
been saved had the latest version of the documentation accompanied the hardware. 
Much time was 6pent deciphering inconsistencies and errors that had already been 
fixed in later versions of the documentation. Also, it was known that there was a 
flaw in the circuit board output buffer. It should be discussed in an errata appendix 
to the circuit board documentation. 

4 Provide applications notes (tutorial) on how the system works and interfacing 
requirements such as timing, reset sequences, andTTL to CMOS voltage interfacing. 
It is likely that many of the people wanting to use these chips will be 
communications oriented, not digital circuit designers. Some documentation should 



be directed to general system design and construction, requirements affecting 
incorporation of these chips into other systems. 

5 Provide documentation on interleaving and deinterleaving, perhaps includ in g 
advantages/disadvantages, implementation techniques and a sample circuit to work 
with this system. 
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APPENDIX 7.A 
Experimental Data 

Test data were designed to subject the syBtexn to a variety of conditions. In 
general the data were designed to provide a large number of code blocks with 
correctable errors, either multiple short bursts scattered in each code block, or a few 
long bursts more widely spaced. The following table shows distributions of some of 
the tests run on the system. 


TABLE 7.1 Test Data Statistical Distributions 


Burst Length 
Mean Variance 
(Gaussian) 

1 15 

2 15 

4 12 

6 12 

8 8 

10 8 

1 7 

2 7 

4 6 

6 6 

8 4 

10 4 


Burst Spacing 
T.nmhH« 

(Poisson) 

255 

265 

255 

255 

255 

255 

128 

128 

128 

128 

128 

128 


The following data show details of tests 4, 9, and 12 as discussed in the results 
section above. Histograms document the distribution of the random errors applied 
to the test apparatus. Note that the graphs are connected point-to-point, not as 
stairsteps. 
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TABLE 7.2 Test 4, GausBian Burst Length. Data 


Burst 

Length 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


N umb er of 
Occurrences 
867 
2043 
4349 
8017 
13143 
19133 
24506 
27689 
27676 
24487 
19131 
13205 
8024 
4347 
2041 
867 
318 
113 
35 
6 
3 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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TABLE 7.3 Test 9, Gaussian Burst Length Data 


Burst Number of 

Length Occurrences 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


39 

223 

979 

3317 

8799 

18387 

29989 

38287 

38309 

29994 

18339 

8808 

3283 

972 

232 

37 

3 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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TABLE 7.4 Test 12, Gaussian Burst Length Data 


Burst N umb er of 

Length Occurrences 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


45637 

45443 

39543 

29736 

19590 

11046 

5492 

2352 

877 

293 

66 

20 

5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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Figure 7.1 Test 4 Gaussian Distribution, Mean = 1, Variance 
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Burst Length 

Figure 7.2 Test 9 Gaussian Distribution, Mean = 8, Variance 
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0.00 37.5 0 75.00 112.50 150.00 197 50 225 00 262.50 

Spacing Between Bursts *10 

Figure 7.6 Poisson Distribution, Lambda = 266 



APPENDIX 7.B 
Program Listings 


This appendix con tains printouts of the programs developed to control an d 
supply data to the encoder/decoder. There were many variations of the immedi a t e 
mode programs, applying various control sequences and single bit or byte errors to 
different parts of the data blocks. For example, the reset sequence in QUT1 9F - A S M 
was put into a loop to debug a problem with initilizing the system. Progressively 
moving a single bit error to new locations on the data block by adjusting loop 
counters permits identification of timing marks on the oscilloscope. 

OUT19F.ASM This is an Immediate mode program used in in i ti al testing. It 
performs the following functions: 

• clears data and error ports 
- resets encoder and decoder 
. then begins infinite loop to output 

— alter nating zeros end ones data, inserting an er ror on the 218 th data byte 

- the check bytes. 

BURST.FOR Creates the Gaussian burst data of a specified mean and variance, 
storing the data in file BURST.DAT. 

STAT.FOR Calculates the sample mean and variance for the data in BURST.DAT. 
Note that the data is constrained, there can be no bursts of length less than zero. 


46 


47 


The upper magnitude limit is 2 16 , since the data is represented as a two byte integer. 
This means that although the user specifies a desired mean and variance in 
BURST.FOR, the experimental values may be slightly skewed. STAT.FOR calculates 
the actual experimental statistics. 

POISSON.FOR Creates the Poisson spacing data given the input X = 1.0/T, where 
T is the spacing between bursts. The data is stored in the file POISSON.DAT 

COUNT2.FOR Creates the loop counter file COUNT.DAT. These are the number 
of bytes of corrupted and uncorrupted data needed to fill out the output array. 

COUNT3.FOR Creates the loop counter files LOOP1.DAT, and LOOP2.DAT. These 
counters are associated with the arrays BURSTO, and POISSON!] that are 
continually filled and emptied of the data being read from the RAM disk. 

CPROG.EXE This program will read through the BURST.DAT, and POISSON .DAT 
data files o nce to calculat e the number of code hlnrlra r nntjining 17 or more errors 
per 265 bytes. 

OUT22.ASM Uses all of the data created above to drive the R/S circuit board with 
Poisson distributed error bursts having Gaussian distributed lengths. 



7.BJ. Program OUT19F.ASM 


port on I/O card ; 

data .port = 0X0378 = 888; 
statui_part = 0X0379 = 889; 
controi_port = 0X037 A = 890; 
port.control = 0x037B = 891; 
port on video card 

data_port s OzOSBC = 956 
statua_port = OxOSBD = 957 
controLport = OzOSBE = 958 
port.control = OrOSBF = 959 


; Static Name Aliases 
TITLE out3.c 
NAME out3 

.286p 

.287 

.TEXT SEGMENT WORD PUBLIC ’CODE’ 

* r ky r 

IdATA SEGMENT WORD PUBLIC DATA’ 

.DATA ENDS 

CONST SEGMENT WORD PUBLIC ’CONST 
CONST ENDS 

_BSS SEGMENT WORD PUBLIC ’BSS’ 

JBSS ENDS 

DGROUP GROUP CONST, .BSS, .DATA 

ASS UME CS: .TEXT, DS: DGROUP, SS: DGROUP 

;EXTRN acrtusedABS 

; EXTRN _chkstfcNEAR 
.TEXT SEGMENT 

ASSUME CS; .TEXT 
; Line 7 

PUBLIC .main 
.main PROC NEAR 
; CLEAR AX REGISTER 
zor az^tz 


; DISABLE INTERUPTS, CLOCK 
cti 

;CLEAR DATA PORT AND ERROR PORT 
mov dz388 
mov ax,Q 

out dx,al ; clear data port 

mov dx4)56 
oat dx,al 

mov dz390 ; clear error port 

mov az4.5 

out dxxl ; clear control portx yarstn high 
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.CLEAR ENCODER BY CLOCKING 32 TIMES WITH INC LOW, LEAVE INC LOW 
;AND CLOCK HIGH ** REMEMBER PORT 890 HAS INVERTED LOGIC** 

mov cx,32 
mov <1x4190 
clkO: 

mov ax,15 

out dbwd ;CLOCK LOW & INC LOW till 

mov ax43 
mov axj.3 

out dx^l ;CLOCK HIGH & INC LOW 1101 

loop cikO 

; RESET DECODER BY CLOCKING 8 TIMES WITH SysRstN LOW, THEN 7 
•.TIMES WITH SysRstN HIGH 

mov ax. 9 BRING SysrstN LOW 1101 

out dx.al 

mov cx,9 
clkl: 

mov ax,ll 
out dx,al 
mov ax,9 
mov ax,9 
mov ax,9 
out dx,al 
loop clkl 

mov ax, 13 

out daual ;BRING SysrstN HIGH 1001 

mov cx,9 
clk2: 

mov ax, 15 
mov 8x4-5 
out dx,al 
mov ax, 13 
mov ax, 13 
mov ax4.3 
out dx^l 
loop dk2 

;CLEAR DATA PORT AND ERROR PORT 
mov dx4388 
mov ax,0 

out dx,al ; clear data port 

mov dx,966 
mov ax,0 
out dx,al 


•.CLOCK 7 TIMES 
;CLOCK LOW 1011 

;CLOCK HIGH 1001 


;CL0CK 8 TIMES 
;CLOCK LOW 1111 

•.CLOCK HIGH 1101 


; clear error port 


mov dx890 
mov axfi 

out dxai ;lnc high & dock high 0001 


start: ; START MAIN LOOP 

; THE FOLLOWING LINES ARE COMMENTED OUT 
mov dx888 
mov ax2£5 

; out dxai ;put ail ones out first; 


mov dx,890 
mov ax, 7 
out dxai 
mov ax5 
out dxai 
mov ax, 7 
out dxai 
mov ax.5 
out dxai 


;inc high & dock low 0011 
;inc high & dock high 0001 
line high Sc dock low 0011 
;mc high Sc dock high 0001 


;OUTFUT 220 DATA BYTES ALTERNATING ZEROS 00000000 AND ONES 11111111 
;WITH INC HIGH 


mov cxllO 
dump: 


mov 

dx888 

mov 

axO 

out 

dxai 

mov 

dx890 

mov 

ax7 

out 

dxai 

mov 

axA 

mov 

axfi 

out 

dx^l 

mov 

dx888 

mov 

ax25fi 

out 

dxai 

mov 

dx890 

mov 

ax7 

out 

dxai 

mov 

axfi 

mov 

axfi 

out 

dxai 

loop dump 


;SET LOOP COUNTER 


;INC HIGH Sc CLOCK LOW 0011 
•.SOMETIMES USED TO PROVIDE DELAY 

;INC HIGH Sc CLOCK HIGH 0001 

;LOAD ALL ONES INTO DATA REGISTER 


jINSEBT ERROR ON THE LAST DATA BYTE 
mov dx£56 
mov ax2&fi 

out dxai -.OUTPUT THE ERROR 11111111 


mov dx£90 
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mov ax,7 

out dxmi ;INC HIGH & CLOCK LOW 0011 

mov 

out dx^I ;INC HIGH & CLOCK HIGH 0001 

mov dx4)56 
mov ax£55 

out dz^al ; output the error 00000000 

mov dz390 
mov ax,7 

out dx^i unc high & clock low 0011 

mov ax,6 

out dx^l 3 dc high & clock high 0001 

mov dx,956 

mov ax^55 

out dx,al 

mov dx.890 

mov ax,7 

out dx.al ;inc high & clock low >3011 

mov ax,5 

out dx,al ;inc high & clock high 0001 

mov dx,956 ; clear errors 

mov ax,0 
out dx^l 

•.OUTPUT CHECK BYTES BY CLOCKING 32 TIMES WITH INC LOW 

mov cx,32 ;SET LOOP COUNTER 

mov dx$9Q 
mov ax4-3 

out dxxl ; bring inc low 

clk3: 

mov ax45 

out dx^l ;INC LOW & CLOCK LOW 1011 
mov ax43 
; mov ax43 

out dx^l ;INC LOW & CLOCK HIGH 1001 

loop clkS 
mov ax£ 

out rix.al ; bring inc high 

; movdx,966 
mov ax,0 
out dxml 
loop start 

.main ENDP 
_TEXT ENDS 
END 


ooooooooooooooo 


7.B.2 Program BURST.FOR 


PROG RAM BURST JOE 

THIS PROGRAM GENERATES AN N SAMPLE GAUSSIAN DISTRIBUTED 
SEQUENCE WITH GIVEN MEAN AND VARIANCE 

INPUT PARAMETERS 

N NUMBER OF REAL WHITE NOISE SAMPLES DESIRED 
VAR DESIRED VARIANCE 
MEAN DESIRED MEAN 

ISEED SEED NUMBER FOR RANDOM NUMBER GENERATOR 
OUTPUT PARAMETERS 
W OUTPUT ARRAY 


REAL W(100), MEAN, VAR 

INTEGER N, ISEED 

OPEN (UNIT =7,FI LE=D:\ BURST JD AT) 

PRINT*, ’ENTER THE NUMBER OF SAMPLES TO GENERATE’ 

READ* N 

PRINT*,” 

PRINT*, ’ENTER THE MEAN’ 

READ*, MEAN 
PRINT*,” 

PRINT*, ’ENTER THE VARIANCE OF THE NOISE S AMPLE S’ 

READ*, VAR 

PRIN T*,” 

PRINT*, 'ENTER THE SEED NUMBER OF RANDOM GENERATION’ 
READ* ISEED 

CALL WGN(N,VARJLSEED,W,MEAN) 

STOP 

END 

C ** ***** i !■**•••••••• WGN SUBROUTINE *<**••••••*********** 

SUBROUTINE WGN(N,VARJSEED,W,MEAN) 

REAL MEAN, W(l), TWOPI, ROOTVA, Ul, U2, TEMP, RAN 
INTEGER I, K, L, M 
INTEGER J 
R=0 

TWOPI=8*ATAN(l .0) 

MslOO 

ROOTVA = SQRTCVAR) 

IF (MODCN42) .NE. 0) M=N+1 
5 CALL RANDOM(MJSEED,W) 

L s 60 
DO 10 1=1 L 
U1=W(2*I-1) 

U2=W(2*D 

IF (Ul XK 0.0) Ul=0.5 


TEMP=SQRT(-2.0*ALOG(U1 )) 

W(2*I-1 )=TEMP*C0S(TW0PI*U2)*R00TVA 
W(2*D=TEMP*SIN(TWOPI*U2)*ROOTVA 
CONTINUE 
DO 30 1=14.00 

W(I)= W (D+MEAN 

IF ((W(D .LT. 256) .AND. (W(D .GE.O)) THEN 
J=INTCWO)) 

K=K+1 

WRITE (7,45) J 
IF (K .EQ. N) GO TO 40 
END IF 
CONTINUE 
IF (K.LT.N) GOTO 5 
CLOSECUNIT=7) 

F0RMATO2) 

PRINT*. ’NUMBER OF SAMPLES GENE RATED’, K 
RETURN 
END 

«*»***»****•*.» RANDOM SUBROUTINE ♦**»*****«*»« 
SUBROUTINE RANDOM(M,ISEED,W) 

REAL W(l), RAN 
DO 20 1=1 M 

W(D=RAN(ISEED) 

CONTINUE 

RETURN 

END 

****«*•»**•«*...**» FUNCTION »«»««••••*««»** 
FUNCTION RANGSEED) 

INTEGER A,C 
A=25173 
C=13849 
M=65536 

ISEED=MOD(A*ISEED+C,M) 

RAN=FLOAT(ISEEDVFLOAT(M) 

RETURN 

END 



oooooooooooooooo 


7.B.3 Program POISSON.FOR 


PROGRAM POISSON.FOR 

THIS PROGRAM GENERATES AN N SAMPLE NORMALLY DISTRIBUTED 
NOISE SEQUENCE AND AN N SAMPLE SEQUENCE OF POISSON 
SPACED NUMBERS 

INPUT PARAMETERS 

N NUMBER OF REAL NOISE SAMPLES DESIRED 
ISEED SEED NUMBER FOR RANDOM NUMBER GENERATOR 
LAMBDA AVERAGE POISSON SPACING 

OUTPUT PARAMETERS 

W REAL ARRAY OF DIMENSION Nxl OF NOISE SAMPLES 
P REAL ARRAY OF POISSON SPACED NUMBERS 


REAL LAMBDA P, RANDOM, W 
INTEGER N, I, A C, M, K, ISEED 
INTEGER J 

PRINT*, 'ENTER THE NUMBER OF SAMPLES TO GENERATE’ 

READ*, N 

PRINT*,” 

PRINT*, TENTER THE AVERAGE POISSON SPACING’ 

READ*, LAMBDA 
PRINT*,” 

PRINT*, TENTER THE SEED NUMBER OF RANDOM GENERATION’ 
READ*, ISEED 

LAMBDA = 1.0/LAMBDA 

A=25173 

C= 13849 

M=66S36 

fcO 

OPEN (UNrr=7,FrLE=T) POISSON.DAT) 

10 ISEED=MOD(A*ISEED+CM) 

W=FLOAT(ISEEDVFLOAT(M) 

P=<-1 .0*ALOG(1 .0-W)VLAMBDA 
IF (P .GE. 0.0) THEN 
J = INTCP) 

K = K + 1 
WRITE (7,20) J 
ENDIF 

IF (KLT.N) GOTO 10 

PRINT*, ’NUMBER OF SAMPLES GENERATE D\K 
20 FORMATCI4) 

STOP 

END 



7.B.4 Program COUNTI.FOR 


PROGRAM COUNT1.FOR 

C THIS PROGRAM CREATES THE COUNT .DAT WHICH CONTAINS 

C THEN NUMBER OF SAMPLES TO READ FROM BURST AND POISSON 
C TO FILL THE OUTPUT ARRAY OF APPROXIMATE SIZE GIVEN 
C 

C ARRAY=255*32=81 60 

C 

C OUTPUT FILE IS COUNT.DAT 

C ALSO THE PROGRAM WILL OUTPUT J, THE NUMBER 

C OF ELEMENTS IN COUNT.DAT 
C 

INTEGER ARRAY, Y,COUNT(20000) 

INTEGER B,P 

PRINT*, ’ENTER THE NUMBER OF SAMPLES TO COUNT 
READ*, N 

PRINT*, ’ENTER THE ARRAY SIZE (OUTPUT) MULTIPLE OF 255 (8160)’ 
READ*, ARRAY 

OPEN (UNIT=5,FILE=D:\BURSTJDAT) 

OPEN (UNIT=7,FILE=D:\POISSONDAT) 

Y=0 

K=0 

J=0 

L=0 

DO 20 1=1 ,N 

READ (5,*) B 
READ (7,*) P 
Y=Y+B+P 
K=K+1 

IF (Y .GE. ARRAY) THEN 
J=J+1 

COUNT J)=K 
Y=0 
K=0 
ENDIF 
20 CONTINUE 

IF (K .NE. 0) THEN 
J=J+1 

COUNT J)=K 
ENDIF 

CLOSE (UNIT=5) 

CLOSE (UNIT=7) 

OPEN (UNIT=5 ,FELE =’C : \ PROG\ MINE \ C OUNT.DAT) 

DO 60 1=14 

WRITE (5,40) COUNT!) 

60 CONTINUE 
40 FORMAT (15) 

PRINT*, ’NUMBER OF ELEMENTS IN COUNT.DAT J 

STOP 

END 


7.B.5 Program COUNT2.FOR 


PROGRAM COUNT2TOR 

C THIS PROG READS FILE COUNT.DAT AND CREATES 
C THE FILES LOOPl.DAT AND LOOP2.DAT 

C 

C ARRAY=3400 

INTEGER B(20000),Y,Z,LOOP1(1000)JjOOP2(1000)ARRAY 
PRINT*, 'ENTER THE NUMBER OF ELEMENTS IN COUNT.DAT 
READ*,N 

PRINT*, ’ENTER THE ARRAY SIZE (BURST & POISSON)’ 

READ*, ARRAY 

OPEN (UNIT=5 ,FILE=’C : \ PROG\ MINE\ COUNT.DAT) 

DO 10 1=1, N 

READ (6,*) B(I) 

10 CONTINUE 

CLOSE (UNIT=5) 

Y=0 

K=0 

z=o 

J=0 

DO 20 1=1 N 
Y=Y+B(I) 

Z=Z+B(I) 

K=K+1 

IF (Y .GE. ARRAY) THEN 
J=J+1 

LOOPl(J)=K 
LOOP2(J)=Y 
Y=0 
K=0 
ENDIF 
20 CONTINUE 

IF (K .NE. 0) THEN 
J=J+1 

LOOPl(J)=K 

LOOP2(J)=Y 

ENDIF 

OPEN (UNIT=7,FILE=’C:\PROG\MINE\LOOPl.DAT) 

DO 30 1=1 J 

WRITE (7,100) LOOPKD 
30 CONTINUE 
CLOSE (UNTT=7) 

OPEN (UNIT=7,FILE=’C:\PROG\MINE\LOOP2.DAT) 

DO 40 1=1 J 

WRITE (7,100) LOOP20) 

40 CONTINUE 

PRINT*, 'NUMBER OF ELEMENTS IN LOOPl AND LOOP2’J 
100 FORMAT (16) 

STOP 

END 


7.B.6 P r ogr am CPROG.FOR 

C THIS PROG CALCULATES THE NUMBER OF UNCORRECTABLE 

C CODE WORDS IN THE DATA GENERATED 
C 

INTEGER POISSO(3500),BURSTC3500),COUNT(6500) 

INTEGER LOOP1(100),LOOP2(100).OUTPUTC8200),TEMP(256) 
INTEGER R£,T,W,X,M.NY,ZJ,K,STOP, START 
PRINT*, 'ENTER THE NUMBER OF ELEMENTS IN COUNT.DAT 
READ*, NCOUNT 

PRINT*, 'ENTER THE NUMBER OF ELEMENTS IN LOOPl.DAT 
READ*, NLOOP 

OPEN (UNIT=5,FILE=’C:\PROG\COUNT.DAT) 

DO 10 1=1 NCOUNT 

READ (5,*) COUNTII) 

10 CONTINUE 
CLOSE (UNIT=5) 

OPEN (UNIT=5.FILE=’C:\PROG\LOOPl.DAT) 

DO 20 I=lvNLOOP 

READ (5,*) LOOPKI) 

20 CONTINUE 

CLOSE (UNIT=5) 

OPEN (UNIT=5,FILE=’C:\PROG\LOOP2JDAT) 

DO 30 1=1 ,NLOOP 

READ (5,*) LOOP2(I) 

30 CONTINUE 

CLOSE (UNIT=5) 

R=0 

T=0 

s=o 

K=1 

NN=31 

OPEN (UNIT=5,FILE=T):\BURSTJ3AT) 

OPEN CUNIT=7,FILE=D:\POISSONJDAT) 

OPEN (UNTT=6,FILE=’ERROR.DAT) 

DO 40 Z=l, NLOOP 
C PRINT*, Z 

W=1 

DO 60 N=1L00P2(Z) 

READ (6,*) BURSTIN) 

60 CONTINUE 

DO 60 N=1L00P2(Z) 

READ (7,*) POISSO(N) 

60 CONTINUE 

DO 70 N=1L00P1(Z) 

Y=1 

DO 80 X=l, COUNT! K) 

1=1 

90 IF (I .LE. BURSTfW)) THEN 

OUTPUT(Y)=255 

C PRINT*, Y,OUTPUT(Y) 


58 


Y=Y+1 

IF (Y .GT. 8160) I=BURSTCW) 

1 = 1+1 
GO TO 90 
END IF 
1=1 

110 IF Q .LE. POISSOCW)) THEN 

OUTPUT(Y)=0 

C PRINT*, Y.OUTFUTCY) 

Y=Y+1 

IF (Y .GT. 8160) I=POISSO(W) 

1 = 1+1 

GO TO 110 
ENDIF 

C PRINT*, BURST(W),POISSO(W) 

W=W+1 

80 CONTINUE 

IF (N .EQ. 1 .AND. Z .EQ. 1) GO TO 135 
S=0 

R=R+1 

TEMP(254)=OUTPUT(l ) 

TEMP(255)=OUTPUT(2) 

DO 130 X=l,255 

IF (TEMPGO .EQ. 255) S=S+1 
130 CONTINUE 

IF (S .GE. 17) T=T+1 
135 DO 140 X=l,253 

TEMP(X)=OUTPUTC7907+X) 

140 CONTINUE 

DO 150 X=l,7905 

OUTPUT(X)=OUTPUT(X+2) 

150 CONTINUE 

IF (Z .EQ. NLOOP .AND. N .EQ. LOOPl(Z)) NN=28 
DO 160 1=1, NN 
R=R+1 
S=0 

START=((M)*255)+1 

STOP=I*255 

DO 170 M=START,STOP 
IF (OUTPUTCM) J5Q. 255) S=S+1 
170 CONTINUE 

IF (S .GE. 17) T=T+1 
160 CONTINUE 

K=K+1 

C PRINT*, N 

70 CONTINUE 

WRITE (6,*) Z,R,T 
40 CONTINUE 
STOP 
END 


7.B.7 Program QUT22.C 


; port on I/O card ; 

; data_port = 0X0378 = 888; 

; statu8_port = 0X0379 = 889; 

; control _port = 0X037A = 890; 

; port .control = 0x037B = 891; 

; port on video card ; 

; data_port = OzOSBC = 956; 

; status_port = OxOSBD = 957; 

; controLport = 0r03BE = 958; 

; port_control = 0x03BF = 959; 

; Static Name Aliases 

i 

TITLE c:\prog\diskl3.c 
NAME diskl3 


.286p 

.287 

_TEXT SEGMENT WORD PUBLIC ’CODE’ 

.TEXT ENDS 

.DATA SEGMENT WORD PUBLIC DATA’ 

.DATA ENDS 

CONST SEGMENT WORD PUBLIC ’CONST 

CONST ENDS 

_BSS SEGMENT WORD PUBLIC BSS’ 

_BSS ENDS 

DGROUP GROUP CONST, _BSS, .DATA 

ASSUME CS: .TEXT, DS: DGROUP, SS: DGROUP 

EXTRN acrtusediABS 

EXTRN _fdose:NEAR 

EXTRN _ent:NEAR 

EXTRN JopemNEAR 

EXTR N _chkstfcNEAR 

EXTRN JscanfcNEAR 

EXTR N .gets :NEAR 

EXTRN .printfcNEAR 


.DATA 

SEGMENT 

$SG223 

DB 

V, 00H 

$SG224 

DB 

’ count, daf, 00H 

$SG228 

DB 

’%u\ 00H 

$SG229 

DB 

Y, 00H 

$SG230 

DB 

'loopl.daf, 00H 

$SG234 

DB 

’%u ’, 00H 

$SG235 

DB 

Y, 00H 

$SG236 

DB 

’loop2.daf, 00H 

$SG240 

DB 

’%u’, OOH 

$SG241 

DB 

Y, OOH 

SSG242 

DB 

’d:burst.dat’, OOH 

$SG243 

DB 

Y, OOH 
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$SG244 

DB 

’d:poisson.dat\ 00H 

SSG248 

DB 

’%u\ 00H 

$SG252 

DB 

’%u\ 00H 

$SG256 

DB 

’%u’, 00H 

.DATA 

ENDS 


.TEXT 

SEGMENT 


ASSUME 

CS: .TEXT 

; Line 10 

PUBLIC 

_main 


_main PROC NEAR 
xor ax, ax 

; DISABLE INTERUPTS, CLOCK 
cli 

•.CLEAR DATA PORT AND ERROR PORT 
mov dx,888 
mov ax,0 

out dx,al ;dear data port 

mov dx,956 
out dx,al 

mov dx,890 ;clear error port 

mov ax4.fi 

out dx,al -.clear control port,sysrstn high 

;CLEAR ENCODER BY CLOCKING 32 TIMES WITH INC LOW, LEAVE INC LOW 
UVND CLOCK HIGH **REMEMBER PORT 890 HAS INVERTED LOGIC** 

mov cx,32 
mov dxB90 

clkO: 

mov ax4S 

out dx,al ;CLOCK LOW & INC LOW 1111 

mov ax43 

out dx^i ;CLOCK HIGH & INC LOW 1101 

loop dkO 

; RESET DECODER BY CLOCKING 8 TIMES WITH SysRstN LOW, THEN 7 
; TIMES WITH SysRstN HIGH 

mov ax.9 ;BRING SysrstN LOW 1101 

out dx,al 

mov cx,9 

clkl: ;CLOCK 8 TIMES 

mov ax41 

out dx,al ;CLOCK LOW 1111 

mov ax,9 
out dx,al 
loop clkl 


;CLOCK HIGH 1101 


mov ax, 13 

out dx,al ;BRING SysrstN HIGH 1001 


mov cx,9 


clk2: 


mov ax,15 
out dx,al 
mov ax»13 
out dx,al 
loop clk2 


;CLOCK 7 TIMES 
•.CLOCK LOW 1011 
;CLOCK HIGH 1001 


•.CLEAR DATA PORT AND ERROR PORT 
mov dx,888 


mov 

out 

mov 

mov 

out 

mov 

mov 

out 


ax,0 
dx.al 
dx,956 
ax.0 
dx,al 
dx,890 
ax, 5 
dx,ai 


;clear data port 


•.clear error port 


;inc high & dock high 0001 


; ARRAY AND VARIABLE INTERNAL SETUP 

push bp 
mov bp,sp 
mov ax ,25658 

call chkstk 

push si 

; poisson = -25418 

; loopl = -18418 

loop2 = -25618 
; w = -8204 

x = -11708 
; m = -25620 

; n = -26634 

; y = -18214 

; z = -18216 

i s -18212 
; k = -18218 

; stop = -8202 

; start = -25632 

; count = -18210 

; zurst = -11706 

"T output = -8200 

; PAUSE = -25630 

infilel = -8206 
; infile2 = -11710 

; Line 19 



push 

push 

call 

add 

mov 

■ ; Line 20 

mov 

SL20001: 

; Line 21 

; Line 22 

mov 

lea 

push 

push 

push 

call 

add 

: Line 23 
inc 
cmp 
jbe 

; Line 24 

push 

call 

add 

; Line 25 

push 

push 

call 

add 

mov 

; T.ina 26 

mov 

$L20002: 

; Line 27 

; Line 28 

mov 

shl 

lea 

push 

push 

push 

call 

add 

; T,in<» 29 
inc 
cmp 
jbe 

; Line 30 

push 

call 


OFFSET DGROUP:$SG223 
OFFSET DGROUF.SSG224 
_fopen 
sp,4 

WORD PTR [bp-8206], ax ;infilel 

WORD PTR [bp-25634], 1 ;n 


si, WORD PTR [bp-256341 ;n 

ax,[bp-18210][si] 

ax 

OFFSET DGROUP:$SG228 
WORD PTR [bp-8206] ;infilel 
_fscanf 
sp,6 

WORD PTR [bp-25634] ;n 

WORD PTR [bp-25634],6277 ;n 

SL20001 

WORD PTR [bp-8206] ;infilel 

_fdose 

sp,2 

OFFSET DGROUP: SSG229 
OFFSET DGROUP:$SG230 
_fopen 
sp,4 

WORD PTR [bp-8206], ax ;infilel 

WORD PTR [bp-26634], 1 ;n 


si,WORD PTR [bp-25634] ;n 

siJL 

ax,[bp-l 841 8][si] 
ax 

OFFSET DGROUP:$SG234 
WORD PTR [bp-8206] ;infilel 
_fscanf 
sp,6 

WORD PTR [bp-25634] ;n 

WORD PTR [bp-25634], 69 ;n 

$L20002 

WORD PTR [bp-8206] ;infilel 

JFdose 
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add sp,2 
; Line 31 


push 

OFFSET DGROUP:$SG235 

push 

OFFSET DGROUP:$SG236 

call 

_fopen 

add 

sp,4 

mov 

WORD PTR [bp-8206], ax ;infilel 

; Line 32 

mov 

WORD PTR [bp-26634], 1 ;n 

$L20003: 
; Line 33 
; Line 34 

mov 

si, WORD PTR [bp-25634] ;n 

shl 

sij. 

lea 

ax,[bp-25618][si] 

push 

ax 

push 

OFFSET DGROUP:$SG240 

push 

WORD PTR [bp-8206] ;infilel 

call 

_fscanf 

add 

sp,6 

; Line 35 

inc 

WORD PTR [bp-26634] ;n 

cmp 

WORD PTR [bp-25634], 59 ;n 

jbe 

$L20003 

; Line 36 

push 

WORD PTR [bp-8206] ;infilel 

call 

_fclose 

add 

sp,2 

; Line 38 

lea 

ax, WORD PTR [bp-25630] ;PAUSE 

push 

ax 

call 

_gets 

add 

sp,2 

; Line 39 

mov 

WORD PTR [bp-18218], 1 ;k 

; Line 40 

push 

OFFSET DGROUP:$SG241 

push 

OFFSET DGROUP:$SG242 

call 

_fopen 

add 

sp,4 

mov 

WORD PTR [bp-8206], ax ;infflel 

; Line 41 

push 

OFFSET DGROUP:$SG243 

push 

OFFSET DGROUP: $SG244 

call 

Jfopen 

add 

sp,4 

mov 

WORD PTR [bp-11710], ax ;inffle2 

; Line 42 

mov 

WORD PTR [bp-1 821 6], 1 ;z 

jmp 

$F245 

SL20004: 
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; Lane 47 
; Line 48 

mav si, WORD PTR [bp-25634] ;n 

lea ax,[bp*11706][si] 
push ax 

push OFFSET DGROUP:$SG252 

push WORD PTR [bp-8206] ;infilel 
call Jscanf 
add sp,6 
; Line 49 

inc WORD PTR [bp-25634] ;n 

$L20005: 

mov si,WORD PTR [bp-18216] ;z 

shl sij. 

mov ax, WORD PTR [bp-25634] ;n 

cmp [bp-256181[si],ax 
jae SL20004 

; Line 50 

mov WORD PTR [bp-2563414- ;n 

jmp SHORT SL20007 

$L20006: 

; Line 51 
; Line 52 

mov si, WORD PTR [bp-25634] ;n 

shl si4 

lea ax,[bp-25418][si] 

push ax 

push OFFSET DGROUP:$SG256 

push WORD PTR [bp-11710] ;infile2 

call Jscanf 

add sp,6 

; Line 63 

inc WORD PTR [bp-25634] ;n 

$L20007: 

mov si,WORD PTR [bp-18216] ;z 

shl si4 

mov ax, WORD PTR [bp-25634] ;n 

cmp [bp-25618][si],ax 

jae $L20006 

; Line 54 

mov WORD PTR [bp-2563414- ;n 

mov si,WORD PTR [bp-18216] ;z 

shl si4 

cmp WORD PTR [bp-1 841 8][si],l 

jae $JCC379 

jmp $FC246 

$JCC379: 

mov si, WORD PTR [bp-18216] ;z 

shl si4 


mov ax,[bp-18418][si] 

mov WORD PTR [bp-25654], ax 



mov si, WORD PTR [bp-8204] 

shi sij. 

lea ax,[bp-25418][si] 
mov WORD PTR [bp-25656], az 
mov ax,WORD PTR [bp-25654] 
mov WORD PTR [bp-25658], ax 
add WORD PTR [bp-25634], ax 
$L20016: 

; Line 65 
; Line 56 

mov WORD PTR [bp-1821 4],1 
; Line 57 

mov WORD PTR [bp-11708], 1 

mov si, WORD PTR [bp-18218] 

mov al,[bp-18210][si] 
sub ah,ah 
cmp ax J. 
jae $JCC449 
jmp SFB262 
SJCC449: 

mov WORD PTR [bp-25644], ax 
mov ax, WORD PTR [bp-26656] 
mov WORD PTR [bp-25646], ax 
mov ax, WORD PTR [bp-25644] 
mov WORD PTR [bp-25648], ax 
add WORD PTR [bp-11708],ax 
SL20010: 

; Line 68 
; Line 69 

mov WORD PTR [bp-18212],l 
; Line 60 

mov si, WORD PTR [bp-8204] 
mov al,[bp-11706][si] 
sub ah, ah 
cmp ax4. 
jb $FB265 

mov WORD PTR [bp-25636], ax 
mov WORD PTR [bp-25638], ax 
$L20008: 

; Line 61 
; Line 62 

mov si, WORD PTR [bp-18214] 
mov BYTE PTR [bp-8200][si],255 
; Line 63 

inc WORD PTR [bp-18214] 

; Line 64 

cmp WORD PTR [bp-18214], 8160 

jbe $1266 
; Line 65 
; Line 66 

mov ax, WORD PTR [bp-26638] 


66 


mov WORD PTR [bp-18212], ax ;i 

; Line 67 
; Line 68 
$1266: 

inc WORD PTR [bp-18212] ;i 

; Line 69 

mov ax, WORD PTR [bp-18212] ;i 

cmp WORD PTR [bp-25636], ax 
jae SL20008 

SFB265: 

; Line 70 

mov WORD PTR [bp-18212], 1 ;i 

: Line 71 

mov bx,WORD PTR [bp-25646] 

cmp WORD PTR [bx],l 

jb SFB269 

mov ax, WORD PTR [bx] 

mov WORD PTR rbp-25640],ax 

mov WORD PTR [bp-25642], ax 

$L20009: 

; Line 72 
; Line 73 

mov si, WORD PTR [bp-18214] ;y 

mov BYTE PTR [bp-8200][si],0 

; Line 74 

inc WORD PTR [bp-18214] ;y 

; Line 75 

cmp WORD PTR [bp-18214], 8160 ;y 

jbe $1270 

; Line 76 
; Line 77 

mov ax, WORD PTR [bp-25642] 

mov WORD PTR [bp-1 8212], ax ;i 

; Line 78 
; Line 79 
$1270: 

inc WORD PTR [bp-18212] ;i 

; Line 80 

mov ax,WORD PTR [bp-18212] ;i 

cmp WORD PTR [bp-25640], ax 

jae $L20009 

$FB269: 

; Line 81 

add WORD PTR [bp-25646], 2 

add WORD PTR [bp-25656], 2 

inc WORD PTR [bp-8204] ;w 

; Line 82 

dec WORD PTR [bp-25648] 

je SJCC636 

jmp $L20010 

SJCC636: 


SFB262: 

; Line 83 


mov 

WORD PTR [bp-18212], 33 

;i 

mov 

WORD PTR [bp-25650], 255 


mov 

WORD PTR [bp-25652], 255 


$L20015: 
; Line 84 
; Line 85 

mov 

ax, WORD PTR [bp-25650] 


sub 

ax,254 


mov 

WORD PTR [bp-25632], ax 

;start 

; Line 86 

mov 

ax, WORD PTR [bp-25650] 


sub 

ax,32 


mov 

WORD PTR [bp-8202], ax 

-.stop 

; Line 87 

mov 

ax, WORD PTR [bp-25632] 

;start 

mov 

WORD PTR [bp-25620], ax 

;m 

jmp 

SHORT SL20012 


nop 

$L20011: 

; Line 88 
; Line 89 

mov 

si, WORD PTR [bp-25620] 

;m 

mov 

al,[bp-8200][si] 


sub 

ah,ah 


mov 

dx,956 


out 

dx, ai 


mov 

dx,890 


mov 

ax,7 


out 

dx,al 


mov 

ax,5 


out 

dx^ai 


; Line 90 

inc 

WORD PTR [bp-25620] 


IL20012: 

mov 

ax, WORD PTR [bp-8202] 

;stop 

cmp 

WORD PTR [bp-26620], ax 

;m 

jbe 

$L20011 


; Line 91 

mov 

ax, WORD PTR [bp-26652] 


sub 

ax,31 


mov 

WORD PTR [bp-26632], ax 

-.start 

; Line 92 

mov 

ax, WORD PTR [bp-25662] 


mov 

WORD PTR [bp-8202], ax 

;stop 

; Line 93 

mov 

ax,WORD PTR [bp-25632] 

-.start 

mov 

WORD PTR [bp-25620],ax 

;m 
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jmp SHORT SL20014 


nop 


$L20013: 
; Line 94 
; Line 95 

mov 

si, WORD PTR [bp-25620] 

;m 

mov 

al,(bp-8200][si] 


sub 

ah,ah 


mov 

dx,956 


out 

dz, al 


mov 

dx,890 


mov 

axj.5 


out 

dx,al 


mov 

axJ.3 


out 

dx,al 


; Line 96 

inc 

WORD PTR [bp-25620] 


SL20014: 

mov 

ax, WORD PTR [bp-8202] 

;stop 

cmp 

WORD PTR [bp-25620], ax 

;m 

jbe 

SL20013 


; Line 97 

add 

WORD PTR [bp-25650],255 


add 

WORD PTR [bp-25652], 255 


cmp 

WORD PTR [bp-25662], 8160 


ja 

$JCC796 


jmp 

$L20015 


SJCC796: 
; Line 98 

inc 

WORD PTR [bp-18218] 

;k 

; Line 99 

dec 

WORD PTR [bp-25658] 


je 

$JCC809 


jmp 

$L20016 


$JCC809: 

$FC246: 

inc 

WORD PTR [bp-18216] 

;z 

$F245: 

cmp 

WORD PTR [bp-18216], 59 

;z 

ja 

$FB247 


; Line 43 
; Line 44 

push 

WORD PTR [bp-18216] 

;z 

push 

OFFSET DGROUP:$SG248 


call 

_printf 


add 

sp,4 


; Line 45 

mov 

WORD PTR [bp-8204], 1 

;w 

; Line 46 



mov WORD PTR [bp-26634], 1 

jmp SL20005 

nop 

$FB247: 

; Line 101 

push WORD PTR [bp-8206] ;infilel 
call _fclose 
add sp,2 
; Line 102 

push WORD PTR [bp-11710] 
call _fclose 
add sp,2 
; Line 103 

push 0 
call _exit 
add ;p,2 
; Line 104 

pop si 
leave 
ret 
nop 

_main ENDP 
_TEXT ENDS 
END 


;infile2 
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l.o ANALYTICAL EXPRESSION FOR THE NEAR HELD 


The near field for a dipole in free space with no surrounding structure is expressed as ( 1): 


E o(r) 


]awUo)_ I z — h e , jkRih + z_+_h_ e _ ]kRvi _ ±L( coskh ) e -jkr 

47 iko \ R\ h Ry, r 


(1) 


and 


Eff\ 


-jwuIAo) 

■ink 


e ~jkR\h 


R \h 


.,-;kRih i 

— — ( coskh ) e~ jkr 

R Vi r 


( 2 ) 


for arbitrary antenna length i ( tip— to— tip) = 2h and k is the complex propagation constant. @ + ja. 

The term fi = — is the phase constant and the term a is the the attenuation constant. The terms 

c 

Rih- ^2h and r are derived from the geometry relationship between the antenna and the observation 
point at which the field is calculated. Figure 1 illustrates the coordinate system and the antenna. 
Cylindrical coordinates are used. The term p (rho) is the radial distance, in the x-v plane, from the 
z axis. The term r is the radial distance from the origin to a point, p = (0. <t>. z) = p(<t>. z) = p(x,y.z). 
Due to antenna pattern symmetry. 6 is not important. 

The electromagnetic Geld calculation code. NEC. will be used to evaluate the possible E and H 
field intensities using a crude model of the Space Shuttle Cargo Bay. To provide a check of the NEC 
code results, an analytical calculation of the E field intensities at one or two points is feasible for the 
dipole in free space. The analytical calculation of the E and H fields for a dipole with the surrounding 
cargo bay support structure is not feasible. However, the estimation of the maximum E field possible 
using a simple “comer reflector ‘ is feasible and will be presented. This will serve to guide 
interpretation of the NEC code. 

As a checkpoint for the NEC code, the E Q and E z field components will be calculated for a free 
space half wave dipole using the analytical expressions, presented as Eq. 1 and Eq. Z for two cases: 
Case I Antenna length t - 7h - A/3 

Case II Antenna length i = 2h = A/2 . 

For each case, the fields will be calculated at an arbitrary point x. y. z of (.6466. 0. -.6466) 
corresponding to the coordinate system illustrated in Figure 1. 

Unfortunately, the text by King and Harrison uses cylindrical coordinates and the NEC code uses 
rectangular coordinates. Fortunately, the conversion between the two is easily accomplished. 


I 
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1.1 Case I Antenna Length, t = lh = A/3 

For this case, it is assumed that the antenna is in free space and hence a = 0. The calculations 
can be accomplished for the dissipative medium of the F layer by specifying the attenuation constant 
a for the F layer. 

The following parameters are assumed: 

Antenna length tip-to-tip. i = 50 meters = 2h 

Frequency, f = 2 MHz: Wavelength. A = j =150 meters 

H = u 0 =4 7T • 10- 7 ; c = 3 • 10 s : k = fi = — = 4.188 - 10‘ 2 . 

c 

using mks units. 

For the point of interest, we have the following values for the parameters indicated: 


o = yX~ + y~ = 6466 

r = + r = v 6466 2 + (- 6466) 2 = .91443 

R\h = 7125 - (- ,6466)| 2 + ,6466 2 = 25.655 
Rvi = J[25 + (- .6466) | 2 + ,6466 2 = 24.362 
/S = — = ~ = 4.189 • 10‘ 2 = .04189 

C A. 

a = 0. 

The one remaining parameter is I^o). the antenna driving point current. The NEC code is run 
with 1 volt applied to the antenna driving point. The resulting current is then calculated by the code 
after calculating the driving point impedance for the specified frequency. 

For the analytical calculation, the driving point current is expressed in King and Harrison but is 
very messy to calculate. Thus, the antenna driving point impedance will be estimated using the ARRL 
Antenna Book curves for a dipole (2) [see Figures 50 and 51 of the next page). 

For the parameters assumed we have as the electrical height in degrees, h = 25 meters = — = 

60 degrees. The resulting estimate for the driving point impedance of a dipole of length l — A/3 
is 


X 

1 


24 - /400 


400.7 


- 86 . 6 ° 


The NEC code has calculated the driving point impedance for the X/3 dipole as (3): 
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Zvec = 26.2 - y427.7 = 428.5 - 86.5 


The two impedances are in fairly good agreement (within approximately 5%). and the NEC code 
impedance should be the more correct as reading the curves from Reference 2 [Figures 49 and 51 
of the previous page] leaves one room for error (the curves are not detailed in fine increments). 

To maintain an independence from NEC code values in the analytical expressions, the estimated 
value for Z will be used. 

The calculated Ir(o) and the NEC code I^o) values are: 

x A 

T i 

l'(o) = — = 2.5 • 10- 3 | 86.6° = 1.48 • KT 4 + y2.495 ■ l(T 3 

esti Z est 

7 — 

3 

and 

l z (o) = 2.33 • KT 3 1 86.49 * = 1.43 • 10" 4 + >233 • 10' 3 

VECi 

T 

Again the results are in good agreement (approximately 6.8%). 

Proceeding to the equation for EJr) , Eq. 1. we have: 

pm _ ;2g2 J 10 6 4jrlQ- 7 2.5 ' 10** 186.6 f - 25.6466 + 24-35 + ^ ecos( 1.097) 

/' 4*4.188- 10' 2 (. 6466) I 25.655 24.36 

3 

or 

Ep(r) = /L16T0 -1 1 86.6* (-.476 + >. 879 + .5231 ->.8517 + .7065 ->.02707) 

A ’ 

= 1.16* 10- 1 1 176.6* [0.7536 -y 2.01- 10" 4 ] = .087417 [ 176.6 • = - .08726 + y.00518 

The £p(r) vector in cylindrical components is equivalent to the following composite sum of 
rectangular coordinates: 

Ejir) = £ t cos <p + E y sirup 


where 4> = tan -1 y/x. 

For our assumed x. y, z coordinates then we have <j> = tan -1 (o) = 0 and 
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Efi) = E, ■ 


Thus, we may compare the calculated EJj) value at the coordinates (.6466. 0. -.6466) directly 
with the Ex value at that point as calculated by the NEC code. 

Eo(?) - - .08726 + y.00518 = .0874171 176.6' 

i' 

T 


E x (r) = -.117 + y.0021 = 0.117| 179.7' 

SEC 

i 

i 

a ■ 034 m 


The result is within a factor of approximately 25% in magnitude. Part ot the contributing factor is 
the use of the actual antenna wire radius in the NEC code, whereas the analytical equation assumes 
an infinitesimally small wire. The NEC code also runs an analysis using the actual optimum antenna 
length versus the theoretical antenna length relative to the wavelength. 

Now let’s calculate the other field component E t (r). From Eq. 2 we have. 


E^ = 


4xk 


e -lkRo, g-iUly, 2 

— — + — (coskh) e~ ,kr 

R\h Rlh r 


-yW4jrlQ- 7 2.5TQ- 3 | 86.6' 
(o/c 


g -i 1.07443 g-yl.02 

25.655 + 


24.362 .91443 


cos( 1.047) e-> 0383 


- y(.075) 1 86.6' (.018563 - ;.0343 + .0215 - y.035 - 1.093 + ;.0419] 
-y(.075) 86.6' (- 1.0529 - j. 0274] = -y(.075) |86.6' ( 1.053256 1 -178.5' 


E£r) = .079 -181.2* = -.0788 + y. 00165 

± L— 

3 

E z (r) = .098 | - 179.29* = - .098 - y.0012 

NEC 

i 

T 

a m 034 m 


Again, we see an agreement to within 19% in magnitude. (The value of a is the antenna wire diameter.] 

The phase associated with the terms is related to the driving point voltage. Reversing the sign 
of the driving point voltage reverses the phase of the NEC voltages. 
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1.2 Case II Antenna Length. l_ =_2A_= A/2 

As a second check, the two field components will be calculated for the half-wave dipole in free 
space at a specific point x. v. z. (.6466. 0. -.6466). 

For this case, we have the following parameters: f = 3 MHz. x. = 100 meters, a = 0: 

„ O) 2jt 71 „ , 

P = — = — = — , h = 25 meters = ‘A • 

c A 50 

The same point for x. y. z will be used and. hence, the values for q, r, Rjh- R-2h will remain the 
same. However, k = 3 is changed. Hence. kRih = 6.2810" 2 (25.655) = 1.61Z kR^h = 1.53 and kr 
= 5.7410- 2 . 

From the ARRL antenna curves of Figures 49 and 50. we have tor this case the estimated driving 
point impedance as 

Z,„ = 2(34.7 - ;20) = “0 + ,40 = S0.6 1 29.7° . 

a/2 

The NEC code calculates 


Z NE c = 80 + ;42.8 = 90.7 28.1 


The agreement is within 11.2% and is considered acceptable. The driving current as calculated 
and as seen by NEC are 


h = 1.24-10- 2 


-29.7* = 1.078 10-2 -;6.1510- 3 


and 


I t = 1.1-10-2 -28.1* 

NEC ^ 

X 

2 


= 9.71- 10- 3 -p.1910- 3 


and agreement is within 12.7% and is acceptable. A 1-volt driving point signal was assumed. 
Using Eq. 1 for E e (r) . we have 


/ 


Ejr) = ;C ^ 0) (- .9997e-^ 612 + .9996e-' 153 4- 1.414cos(1.57)e-> 5 - 74 10 “ ! ) 

= jA6AIjHp) ( + 4.118- 10’ 2 + 7 . 9988 + 4.076- 10" 2 - 7 . 99876 + 1.124- 10~ 3 - y6.459-KT 5 ) 
= .575 | 60.3* (8.310- 2 -y2.46-10- 5 ) = .575 1 60.3* (8.3- 10 ~ 2 [■ITT 1 '! 

E&) = 4.77- 10 “ 2 1 603* = 2.36 10 ~ 2 + >414 10 " 2 

A ' 1 


and 


<1 


E x ( r ) 


vec 

« ui4 m 


9.88- IQ' 2 



-9.2810- 2 + y3.3810' 2 


where a = wire diameter. 
For E t we have 


Eff) = -j3QUo) 

A 


g-; 1.6 12 e -i 1.53 

25.655 + 24.36 


^ ^ ^ -ji.1410- 1 

.91443 


= 3.72- 10- 1 I-119.7 * [- 1.6- 10- 3 + /UfrlO" 2 -t- 1.67- 10" 3 -y4.110- 2 - 1.7510" 3 + jl lQr 4 ] 
= .372 1 - 120- [-1.68-10- 3 -y2.53-10- 2 ] - .372 |-120‘ (2.53- 1C 2 1-93.8 *) 


Erfr) = 9.41- 10~ 3 | - 213.8* = -7.810' 3 + y52310 _3 


and 


Ei 


NEC 

X 


014 m 


115' 10-* | 175.8- = - 114.7- IQ" 3 + ;8.4-10- 3 


1.3 Comments. on Results of Analysis Versus NF.C Code 
It should be noted that the differences between 

E„ and E x and again for E, and E. 

i NEC - SEC 

'■ k 2 i L 
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may be attributed to the .034 m wire diameter used in the NEC code versus the assumed zero diameter 
wire used in the analytical equations. To demonstrate the wire diameter effect, the NEC code was 
rim for various wire diameters. Table 1 compares the results, and the asymptotic behavior of the NEC 
code is illustrated, as is the clear trend towards the mathematical analysis results for zero diameter 
wire. 

To verify the NEC model, the analysis was calculated for A/3 and A/2 lengths as noted previously. 
The NEC code estimates are more realistic and accurate, uti lise numerical solution of integral 
formulations that have no closed form analytic solutions. To check the scaling of results run for 1-volt 
driving point sources, the A/3 free space dipole case was run for 10 volts as well. Table 2 illustrates 
the results. As expected, the E field scales in direct proportion to differences in Ed nvin g po mt and in 
proportion to the square root of the difference in radiated power. / 

This section has presented analytical results and NEC code calculations for companson. At this 
point, one should have confidence that the NEC code results are will be close to those results predicted 
by analysis. The crowning feature of the NEC code is that it may be used to calculate fields for 
structures for which no analytical results are possible, such as the Space Shutde Cargo Bay. 


2.Q NEC CODE RESULTS FOR DIPOT F AND SHUTTLE CARGO 

The field strengths that could be induced in the Shuttle Cargo Bay by the WISP dipole antenna 
are the items of interest. It is not feasible to derive an analytical solution for this structure. However. 
NEC can be used to calculate approximate field intensities in the cargo bay. The term approximate 
is used since it is not feasible to exactly model the Shuttle Cargo Bay in the NEC code. 

An approximate model of the cargo bay may be comprised using a series of “patch plates” and 
an outside closed surface of large rectangular plates as is illustrated in Figure 2a. In Figure 2a. the 
WISP dipole is illustrated as a horizontal dipole lying along the x The Shutde Cargo Bay is 
depicted as lying along the y axis. It should be noted that the orientation of the antenna is in a different 
direction than the orientation used for the c alculati ons of Section L0. This is for convenience in 
modeling the cargo bay with the NEC code. 

From Figure 2a, one may observe the patch structure used to model the cargo bay. The coordinate 
points labeled in Figure 2a are always depicted x, y, z coordinates. It may be noted that the origin 
of the coordinate system is located 21” under the center of the WISP dipole as illustrated in the lower 
right-hand comer of Figure 2a. Figure 2b illustrates the interior of the cargo bay model in greater 
detail. For the NEC code calculation, a series of points in various xjl planes were calculated. These 
points were arbitrarily chosen. The lower figure in Figure 2b depicts the ten calculation points as 
located in the x~z plane for y = 0. The tail of the Shuttle is in the negati ve y direction, and the nose 
of the Shuttle is in the positive y direction. The ten ralmiaM points in the x.z plane were also 
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ILLUSTRATION OF WIRE DIAMETER EFFECT FOR NEC CODE 
ESTIMATE VERSUS IDEAL ANALYTIC SOLUTION 
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As the wire diameter approaches 0, the NEC code incurs numerical limitations. However, the trend towards 
the analytical results is illustrated. 



TABLE 2 

A CHECK ON SCALING 1 VOLT TO 10 VOIT SOURCES 

Antenna is W3 tip-to-tip, f = 2 MHz, X = 150 meters, t = 2h = 50 meters, a = ,034 m 
E fields observed at coordinates *, y, z = (.6466, 0, -.6466). Antenna as Figure I. 

NEC estimated driving point impedance for all cases: z = 26 - j416. 



As expected, the fields scale in direct proportion to Edging point and in proportion to the square root of the 
power difference. 





+ 13.192 


ORIGINAL PAGE IS 

OF POOR QUALITY 





calculated in planes parallel to the x.z plane but for various values of y separated at 1.56 meter 
increments, such as y = -6.0% m. y = -4.536 m. etc., as illustrated in Figure 2b. 

The NEC code results for the Shuttle Loading case are not reliable nor accurate. That this is the 
case is illustrated by the three sets of data in Table 3. Several potential reasons for the nonreliable 
results might be: 

1. Patch size outside window of allowable parameters relative to the wavelength. 

2. Patch coordinates incorrectly specified. 

3. NEC not suitable for such close-in geometry. 

The NEC code problems will be investigated in the following time period of October 1991. 
Hopefully, consistent results can be obtained. 


3.0 CONCLUSIONS 

From Table 2 one may observe that for the HALFWAVE DIPOLE in free space at a point within 
.646 meters of the antenna fairly large field strengths are to be expected. For 2000 volt driving point 
potential 300 watts radiated power, a field strength of 250 volts/m could easily be expected. In the 
presence of the Shutde. the structure could create an increase of 2 or more. Thus, field strengths of 
500 v/m or more could exist in the Shuttle Cargo Bay near the antenna. 
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TABLE 3 


NEC CODE RESULTS FOR SHUTTLE LOADING E-FIELD VALUES 
CALCULATED AT THE POINTS A, B, C. D, E OF FIGURE 2b FOR 

THREE PATCH SIZES 

LACK OF SYMMETRY | E| AT B NOT EQUAL TO | E| AT D, 
AND VARYING RESULTS AT THE SAME POINTS FOR 
VARYING PATCH SIZE INDICATE NON VALID RESULTS , 
^driving point = volts. l A/2. 


PATCH SIZE 





1 1 ■inn 

A 

20 V/m 

18.4 V/m / 

115 V/m 

B 

18 V/m 

13.6 V/m 

10.8 V/m 

C 

15 V/m 

8.6 V/m 

11.20 V/m 

D 

15 V/m 

5.9 V/m 

20.% V/m 

E 

9 V/m 

6.1 V/m 

30.25 V/m 
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